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Abstract 

We investigate the form and stability of a cloud of atoms confined in a 

harmonic trap when the scattering length is negative. We find that, besides 

the known low density metastable solution, a new branch of Bose condensate 

appears at higher density when non locality effects in the attractive part are 

taken into account. The transition between the two classes of solutions as a 

function of the number N of atoms can be either sharp or smooth according 

to the strength and range of the attractive interaction. Use of tight traps is 

favorable for investigating the evolution of the system as the strength of the 

effective interaction increases with N. 
03.75.Fi, 05.30.Jp, 32.80.Pj 
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Few alkali isotopes are characterized by a negative scattering length ar at low energy 
in the triplet channel This is believed to have important consequences on the existence 
and the structure of the Bose Einstein condensed (BEC) state for trapped clouds of atoms 
when the isotope follows Bose statistics. The most celebrated case is 7 Li for which there 
is experimental evidence in favor of a BEC state if the number of atoms does not exceed 
a threshold ||. This behavior was in fact predicted theoretically. When the interaction 
between atoms is represented by an effective local potential v e g(r) = (Auk 2 /m)aT<5 3 (r), a 
uniform gas phase is unstable towards the formation of collapsed state of high density. The 
situation changes in presence of an external trap ||. Within the Gross-Pitaevskii (GP) 
equation for the condensate wavefunction the kinetic energy due to the confinement in 
the trap opposes the collapse as long as the number N of atoms does not exceed a critical 
number N c which depends on the parameters of the system (isotope mass m, scattering 
length ar, shape of the trap, etc.) j5|. Even for N < N c , we must recall that a dilute cloud 
represents a metastable state which can be experimentally observed only if its lifetime is 
long enough. 

The possibility of having positive as well as negative scattering lengths in different atomic 
systems arises because the potential between a pair of atoms generally supports bound states 
leading to the formation of a dimer with several vibrational levels. Under this condition, 
| or | can be very large but at the same time the energy dependent scattering cross section 
(Tt(E) starts to deviate from its zero energy limit ctt(O) = Aira 2 ^ already at an energy that 
can be extremely small 0. The calculations by Cote et al. [|7| for 7 Li show that <jt{E) 
is reduced to 50% of or(0) at energies of the order of 3 10~ 9 a.u. . This value is small 
enough that s-wave scattering is dominant over higher angular momentum states meaning 
that the scattering process between two atoms can still be represented by an effective two 
body interaction i> e g. However, in general we cannot neglect the momentum dependence 
of Ves(k) for the colliding atoms. This amounts to say that the effective interaction is non 
local: In this note we investigate the effects of non-locality on the stability of the cloud. 

It is rather easy to convince that non-locality in the interatomic effective potential can 
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modify the stability condition of the Bose condensate. The GP functional £ with non 
local interaction v c q{y) reads: 

Sm = f d 3 r^\W\ 2 + f d 3 rf/ ext (r) |M/(r)| 2 + 

lJd 3 rJ rfV|vI/(r')| 2 t; cff (r-r')|^(r)| 2 (1) 

where ^(r) is the wavefunction of the condensate and U ext (r) is the potential of the trap. 
In this paper we consider a symmetric harmonic trap U ext (r) = ^moJ^r 2 . In the local limit, 
one recovers the standard form of the GP functional. The ground state wavefunction of a 
cloud of N atoms is determined by minimizing £[*&] with the constraint that 

d 3 r|^(r)| 2 = N (2) 

In the ground state \l/(r) is positive definite and spherically symmetric. 

As a first step, we discuss an approximate, variational approach to this problem which 
already shows the main features of the exact solution. As a trial wavefunction we choose a 
Gaussian form: 
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with a single variational parameter A which defines the size of the cloud in units of the 
harmonic oscillator length an = [h/(muo)] 1 /' 2 . By substituting r' = r + s in Eq. (|l|) and 
developing \l/(r + s) in powers of s, we perform a gradient expansion of the energy functional 
£[$?]■ To lowest order in the gradient of the wavefunction, the energy reads: 

Ar /3 2 3 2 3 5 



5(A) = -^N (J A 2 + t^A + 7 A 3 + 7 T ^ J (4) 

where the two parameters 7 and r depend on the choice of the effective interaction providing 
a dimensionless estimate of the strength 7 and range r of v e g(r): 



7 = J 2/ 7r N ax cljj 

1 -1 



3 / drr v e g(r) 



2a H / drr v eS (r) 



(5) 



According to these definitions, r is basically minus the square of the range of the effective 
interaction and is therefore a negative quantity. For a local interaction r = and we recover 
Baym and Pethick result ||. In this limit and for negative scattering length, 7 < and 
£ (A) goes to —00 for A — > 00 leading to a collapsed ground state. However, S (A) has a local 
minimum for a finite value of A when I7I < 4/5 5//4 , i.e. when N < N c ~ 0.67 a B | | ~ 1 . This 
represents the metastable state mentioned above. 

It is clear from Eq. (§) that non-locality changes the stability condition because, for nega- 
tive scattering length, the product 7T is positive. Depending on the value of the parameters, 
S(X) can have one or two minima and the collapse is prevented in any case. However, this 
result is only suggestive because the gradient expansion breaks down as soon as non-locality 
becomes important and a more appropriate formalism is needed. We have then studied 
the variational equation obtained with the same choice of a Gaussian wavefunction (Q) but 
without invoking gradient expansion. This problem requires the explicit definition of the 
effective interaction. We assume that the attractive potential has a finite range r e and in 
addition we allow for the presence of a repulsive contribution which is modeled as a local 
positive term defined by a "high energy" scattering length a R > 0. With this choice, v e fi(k) 
changes sign from negative at small k (i.e. low energy), to positive at larger k thereby 
mimicking the microscopic computations of the I = phase shift 5k which show a change 
in sign from positive to negative when k increases ||. At still higher momenta, scattering 
in other channels with i 7^ becomes important but we assume that the range of density is 
such that this regime is never reached. 

In the model we study, the effective interaction is then written in the following form: 

Ves(k) = [or + (a T - a R ) f(kr e )] (6) 

m 

We have considered two choices for the shape function f(x): A Lorenzian f(x) = (1 + x 2 )~ 1 
and a Gaussian f(x) = exp(— x 2 ). The results do not depend on the specific choice of f(x) 
and so we will discuss only the Lorenzian case. We use interaction parameters appropriate 
for 7 Li: = —27 a b ||, r e = 10 3 a B and a R = 6.6 a B || (where a B is the Bohr radius). 
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The energy functional for the Lorenzian potential can be analytically expressed in terms of 
elementary functions and reads: 

+r 2 erfc (xX' 1 ) exp ( X 2 A~ 2 ) } (7) 



where 



v , 2 a R 2 a T - a R 

1r= n \\-— Tx = N\\-a H — — 2 — 



7T an V 71" r , 



O'H Ar 2 a T — a R 

X = —f=r- r 2 = Na H 



V2r e " rl 

and erfc(x) = 1 — erf(x) is the complementary error function. The extrema of £(X) are 
obtained as solutions of an algebraic equation. This equation has either one or three positive 
roots depending on the parameters and on number N of atoms in the cluster. When three 
solutions are present, the intermediate one represents an unstable state (i.e. a local maximum 
of the energy) while the other two respectively describe a low density metastable solution 
and a minimum which represents the stable solution within GP approximation. 

The variational results for three typical trap sizes are shown in Fig. 1 where the average 
radius and the density at the center of the cloud are plotted as a function of N. In the same 
figure, the variational data are also compared with the exact solution of the GP equation, 
obtained by numerical integration of the corresponding self-consistent Schrodinger equation. 
In fact, the exact minimization of the GP functional (|1|) gives rise to a non linear eigenvalue 
problem for the ground state wavefunction. The ground state is spherically symmetric and 
can be written as \P(r) = 0(r)/r where 0(r) satisfies: 

d 2 TTIU!^ 

"r) + M0(r) + W{r)(f)(r) = af>(r) (8) 



2m dr 2 2 
Here e is the chemical potential and 



W(r) = J dV^rOlVKr-O (9) 

is the self-consistent interaction. For effective potentials written as sums of Yukawa func- 
tions, this integro-differential equation can be reduced to a boundary value problem for a 



set of ordinary differential equations. In particular, in the case of interest, the form (||) 
gives rise to three ordinary differential equations which have been integrated by standard 
Runge-Kutta method. 

Fig. 1 shows the good performance of the variational approach which is always very 
close to the exact solution. For large N the size of the cloud is remarkably independent of 
the trap size demonstrating that the atoms are in a self-bound configuration. The size of 
the cloud is governed by the effective range: For large values of r e (r e 3> \ar\) an d N not 
too large, the reduced density p|ot| 3 is small. This indicates that the self bound state that 
we find represents a novel regime of the cloud, intermediate between the very low density 
state already predicted within the approximation of local interaction and the collapsed high 
density state which depends on the detailed shape of the true interatomic potential flCfl . 

The asymptotic large N behavior of the exact solution can be obtained analytically from 
our equation: The radius of the cloud reaches a finite limit which coincides with the result 



of the Thomas- Fermi approximation to the GP equation |[L1|| . On the other hand, the low 
density branch is virtually indistinguishable from the exact solution of a local attractive 
potential, also shown in figure. The effects of non-locality become important just when the 
radius of the cloud rapidly drops. This "transition" is discontinuous for large traps, where 
the reentrant behavior of the curve shows the presence of an unstable branch. By reducing 
the trap size, however, this discontinuity is strongly reduced and, below 0.3 /zm, the unstable 
branch disappears and there is a smooth evolution from a very dilute cloud to a less dilute 
state with an increasing density as iV grows. 

The predictions of the standard treatment of clouds of alkali atoms in terms of a local 
pseudopotential have to be modified in two respects when the scattering length is negative. 
In the case of shallow traps, the stability threshold N c is only an upper bound. There is a 
lower threshold N c i for the higher density state and in the intermediate region N c i < N < N c 
the low density branch is metastable towards the higher density solution. For instance, in a 
3 trap N c = 1300 and N c i = 160. Which state is reached experimentally depends on how 
the cloud is formed. In the case of a tight trap a threshold does not exist altogether. This 

6 



feature allows to explore experimentally how the state of the system evolves as the interaction 
strength increases and depletion effects start to set up even at the lowest temperature. It 
should be noticed that the presence of a repulsive short range interaction in our model 
potential does not play a crucial role. The higher density branch as well as the existence 
of two regimes depending on the size of the trap are present even if we put clr = 0. The 
only difference is that, in this case, the radius of the cloud slowly decreases as N — > oo and 
does not reach a finite limit as in Fig. 1. In any case, it should be kept in mind that in the 
large N limit our results are only qualitative because the GP equation itself breaks down 
and interaction effects are expected to produce a depletion of the condensate when p|cit| 3 
is not very small. However, we believe that the trend shown by this equation is significant 
even if quantitative results must wait for a more accurate determination of the interatomic 
potential and a better treatment of of the many body problem. 

Another interesting problem which can be addressed by this formalism is the shape of the 
condensate wavefunction. This information may be useful for the quantitative determination 
of the number of atoms in the condensate from experimental data. The GP ground state 
wavefunction, normalized to unity at the origin, is shown for different numbers of atoms in 
the cloud in Fig. 2a. We see that some deviation from a Gaussian form develops as density 
grows. In particular, a sharp peak at the cloud center grows at large N. Along the very low 
density branch we verified that the wavefunction is not affected by the non locality of the 
interaction and is very close to the numerical result in the local approximation. However, 
significant deviations appear as soon as we approach the edge of the low density branch. 

As a next step we can investigate the dependence of the shape of the equilibrium curve 
when the range of the attractive interaction is varied. In fact, we expect that the larger 
is this lengthscale, the stronger is the effect of non locality in the stability diagram of the 
cloud. More importantly, the value of the effective range r e is rather uncertain even in 
Lithium and it is useful to test the sensitivity of our results to variations of the shape of the 
effective interaction. In Fig. 2b we show the equilibrium curves of a cloud of Lithium atoms 
in a harmonic trap of an = 1 /im when the effective range r e is artificially halved or doubled 
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with respect to the reference value r e = 10 3 ae previously adopted. The results clearly show 
a smoothing of the transition between the two branches induced by the enhancement of non 
locality. By increasing the effective range, the maximum number of atoms in the low density 
branch N c slightly increases and the two branches of solutions tend to merge into a unique 
phase smoothly connecting the low and the higher density regimes. A quantitative analysis 
of the variational solution shows that the occurrence of the single branch regime is mainly 
determined by the ratio between the trap radius and the effective range, independently of the 
scattering length: an < 5r e . Conversely, the number of atoms in the cloud at the inflection 
point can be estimated as N c ~ 4.5 r e |ar| _1 , independently of the trap radius. These results 
show that the effective range plays a crucial role in determining the stability of the bosonic 
cloud at intermediate density. The previous estimates will become particularly valuable 
when an accurate determination of the interaction parameters will become available also for 
other alkali isotopes with negative scattering length. Interestingly, new techniques are being 
devised in order to modify the effective interaction in these systems ||12|| . 

In conclusion, the analysis of the Bose-Einstein condensed ground state of Lithium atoms 
in a harmonic trap has shown that atomic clouds may exist in a novel state, intermediate in 
density between the very dilute regime which has been studied up to now and the collapsed 
high density state which is determined by the short range behavior of the true interatomic 
potential. Also this intermediate state is self-bound because the size of the cloud is essen- 
tially independent of the size of the trap. The importance of the study of this intermediate 
regime arises from the possibility of tuning the strength of the interaction by changing N: In 
this way it will be possible to study how the depletion of the condensate affects the structure 
and the physical properties of the state. To this end, it is important to confine the system 
in sufficiently small traps in order to avoid the discontinuous jump between the two phases. 
Interestingly, such micro-traps are under development [|HJ. Non locality in the effective 
interaction is a crucial ingredient which prevents collapse when the scattering length is neg- 
ative. Long range potentials favor the formation of big clouds which still maintain a rather 
low density. Therefore, it is quite interesting to select elements and isotopes characterized 
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by negative scattering length together with long range effective interactions. We stress that 
our results apply to the dilute regimes where GP equation correctly represents the behavior 
of the ground state. At higher density, collisions between atoms deplete the condensate. 
In addition three body scattering starts to play a role and spin-flip processes lead to rapid 
destruction of the cloud. 
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Figure captions 

Fig. 1 Average radius of the cloud R (panels a,c,e) and central density (panels b,d,f) as a 
function of number of atoms N for the trap sizes an — 3 /im, 1 /im and 0.3 /im respectively. 
Full dots represent the numerical solution of the GP equation with non local Yukawa poten- 
tial. Open dots correspond to the local approximation. Solid line is the variational result 
and the dashed line is the asymptotic limiting radius for N — > oo. Parameters have been 
chosen as to mimic 7 Li. 

Fig. 2 Panel (a): Condensate wavefunction \I/(r) for a trap size a H = 0.3 /im and non 
local interaction. The non interacting result (Gaussian wavefunction) would be a straight 
line in this plot. Dashed line: N = 68 (low density). Dotted line: N = 159 (near the 
maximum of the low density branch). Solid line: N = 166 (higher density branch). Panel 
(b): Equilibrium curves for parameters appropriate for Lithium atoms in a trap of an = 1 /mi 
and three choices of the effective range of the Lorenzian attractive potential. Variational 
results. Dotted line: r e = 500 a B . Solid line r e = 1000 a B - Dashed line r e = 2000 a B - 
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